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I. INTRODUCTION 



The last few years have seen the concept of control landscapes used in connection with a 
variety of problems in quantum control theory [ ]. Generally, a control landscape is a map, 
defined implicitly through the control dynamical system, that assigns a real observable 
value to each admissible control field. In the case of quantum systems, the landscapes 
may be conveniently thought of as the composition of two maps: a control— )-state map and 
a state— j-observable map. The "state" in this case is frequently the unitary time-evolution 
operator U (T, 0) which is the general solution to the Schrodinger equation at some final time 
T; so the control— J-state map takes in a control field and a final time T and produces U (T, 0). 
The state— i-observable map (often called the "kinematic landscape" ) is a real- valued function 
on the special unitary group SU(A^), i.e. it takes in a unitary operator U and produces the 
value of the final observable. Under the composition control— J-state—^-observable (the so- 
called "dynamical landscape"), the goal of optimal control is generally to maximize the 
observable with respect to the control. 

Such landscape formulations have been used to gain some understanding of the critical 
point structure - and thus the nature of the gradient flow - of these optimal control problems 
for state-to-state transitions [2-1], general quantum mechanical observables on an ensemble 
[5-7], and unitary transformation (quantum gate) preparation [8-10]. In addition, these 
studies have considered the implications of the vast multiplicity of controls that are all 
capable of maximizing the landscape value. 

In the present paper, we investigate one aspect of this multiplicity of solutions: the 
topologies of various important subsets of control space. The main result presented here 
is a homotopy equivalence, under certain assumptions, between the set of controls carrying 
one point of a Riemannian manifold to another and the loop space of the reachable set. A 
simple corollary of this result in the case of quantum control is that any connected subset 
of SU(A^) (such as the maximal submanifold of any of the above kinematic landscapes) has 
a connected preimage in control space. The full implications of the wealth of topological 
information offered by the homotopy equivalence are yet to be understood, but this corollary 
on connectedness is important in understanding the gross structure of control landscapes as 
well as the behavior of many optimization and exploration algorithms, especially continuous 
methods such as D-MORPH [11-15]. 

This work has a relation to that of Colonius et. al. [16-18], however it is much more 
closely aligned with the work of Sarychev [19, 20]. Indeed, the arguments presented in this 
paper follow very closely those of Sarychev. Sarychev considered drift-free control problems 
on general Riemannian manifolds and proved a homotopy equivalence between the space 
of trajectories with the usual compact-open (i.e., uniform) topology and the loop space of 
the underlying Riemannian manifold. We extend his results by incorporating a drift term 
and pulling the problem back to the control space with a stronger topology. While we have 
presented this work in terms of quantum control, the principal theorem is proved for affine 
control problems with drift on general Riemannian manifolds. 
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To help place the present work in its proper context, it should be noted that the statement 
of the main theorem in Sarychev's second paper [ ] on the subject is incorrect. As stated, 
the theorem declares that the homotopy equivalence indicated above holds for any "conic 
polysystem" satisfying the Lie algebra rank condition. If true, that theorem would largely 
contain the results of Theorem 1 of the present work. However, Sarychev's proof in [20] 
only covers symmetric conic polysystems, which are systems without drift, suggesting that 
the wording of the theorem was a simple mis-statement. Moreover, Montgomery [21, 22] 
has constructed a counterexample on SO (3) based on Little's [23] work on nondegenerate 
curves on 5*^ which demonstrates that the conclusions of Sarychev's theorem do not hold for 
all conic polysystems. However, Montgomery's example restricts the controls to be strictly 
positive functions. The importance of this difference in the space of controls is made clear 
(as in Montgomery's paper ["'""]) by comparing Little's theorem to Smale's earlier results 
on regular curves [J ;]. Translated into a scalar affine control problem with drift on S0(3), 
Smale's theorem applies to the case in which the control is an arbitrary continuous function 
of time and yields different topology to that of Little's theorem, which applies to the case 
where the control function is constrained to be positive valued. In extending Sarychev's 
work, we consider a class of control problems with drift (thus generally outside the domain 
of the proofs offered by Sarychev) in which the space of controls is taken to be large enough 
to include all piecewise constant functions (thus excluding Montgomery's counterexample). 

The paper is organized as follows. Section II states and proves the main theorem of 
the paper: the homotopy equivalence between the fiber of the endpoint map and the loop 
space of the state manifold. In section III, this theorem is applied to various formulations 
of quantum control to determine the topologies of the corresponding fibers. Methods of 
algebraic topology are employed in Section IV to bootstrap from the topologies of the fibers 
to the topologies of more general preimages of subsets of the state space. The results are 
summarized in Section V. Five appendices are included that prove the continuity of various 
important maps, and describe some similar results that may be obtained for scalar control 
problems on 2-level systems (i.e. qubits). 



II. DYNAMIC HOMOTOPY 



Consider an affine control problem on an n dimensional, second countable, smooth, con- 
nected Riemannian manifold A^, where the control system is of the form 

dz ™ 

- = f,{z) + Y,mm z{Q) = z,eM (1) 

i=l 

and we assume the /j's are smooth, complete vector fields on }A and that the Lie algebra 
generated by the control vector fields {/i, . . . , fm} (excluding the drift /o) spans the tangent 
space of M. at every point. This latter assumption necessarily requires that m be at least 2. 
These assumptions imply that the system is strongly controllable [25, Ch. 4, Thm. 2], i.e. 
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A — * 

that for any pair of points zq, Zf E Ai and any T > 0, there exists a control £ that will steer 
the system from Zq to Zf in some time < T < T. In laboratory implementations of quantum 



control, the m = 1 case of scalar control is more common at present, but experiments with 
m > 2 can be implemented (see also Appendix E). 

Now, let Ko denote the space of (m + l)-tuples (T, £i, . . . , £rn) consisting of a final time 
T > and m functions in L^[0, oo) such that £i{t) = ■ ■ ■ = S^it) = for all t > T. This 
space will be topologized by the metric 



d((Ti,£i,...,£J,(r2,J-i,...,J-J) = |Ti-T2| + ^ / \S,{t)-J^j{t)\dt. (2) 



This topology is somewhat weaker than the corresponding topology, but will be 
necessary for continuity of a certain map later on. We can also define a concatenation 
operation on Kq, denoted by -k, so {Ti,£i,.. . , £m) * (T2, J^i, . . . , Tm) = {Ti + T2, £^1, . . . , ^2), 
where 



Let K denote a subset of Kq that contains the piecewise constant controls, is closed 
under concatenation, and is such that for every element of K, the system (1) admits a unique 
trajectory over the corresponding interval [0, T] for any initial point Zq E M. (any such subset 
will suffice). We will assume also that if C = (T, £^1, . . . , £m) G K, then for any < T' < T, 
the truncation of C to the interval [0,T'] is also an element of K. Now, for any x G A^, 
define a map e^; : K — > in the following way. For a control C = {T,£i, . . . ,£m) G K, 
integrate the control system (1) out to time T with initial condition zq = x and define 
^x{C) := z{T), the solution at the final time T for control C. This will be called the 
endpoint map on the control space. It is a continuous map, as shown in Appendix A. In 
addition e^.(Ci ^ir C2) = c,^(Ci){C2)- 

As indicated in the introduction, the principal result in this paper is the following theorem 
in which PAi is the path space on the state manifold Ai: 

Theorem 1. For each x,y E Ai, the subset t^^{y) G K. of controls carrying x to y is 
homotopy equivalent to QAi, the loop space of the manifold Ai. This homotopy equivalence 
is carried by the restriction of the trajectory map r : K — t- PB, defined implicitly through 
the control system (1) with initial condition x. 

For the sake of clarity, fix now some basepoint b E Ai. The proof of this theorem will 
involve first showing that the endpoint map is a (Hurewicz) fibration and then exploit- 
ing the resulting fibration sequence and corresponding long exact sequence of base-point 
preserving homotopy spaces from some arbitrary pointed topological space. 





(3) 
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A. The Endpoint Map is a Fibration 



We begin with the proof that e?, is a fibration [26, 27]. For a continuous map p : E ^ B 
between topological spaces, let Xp C E x PB consist of all pairs (e, 7) such that p(e) = 7(0), 
where PB is the path space on B with compact-open topology, and let p : PE — )• Xp be 
the map p{ri) = {ri{0),pri). This map p is called a (Hurewicz) fibration if there exists a 
continuous lifting function A : Xp — )■ PE such that po A is the identity on Xp. A continuous 
map p : E ^ B is a local fibration if for each b E B there exists an open neighborhood G 
such that p restricted to p~^{G) is a fibration. The Hurewicz uniformization theorem [26] 
states that if p is a local fibration and B is paracompact, then p is a (global) fibration. Since 
the endpoint map C;, : K — )■ maps to a manifold, the base space is paracompact, so we 
need only prove that Cb is a local fibration to have the desired conclusion. To that end, we 
first prove the existence of a particular cross-section map on a neighborhood of Xq € A^. 

Lemma 1 (The Cross-Section Map). For any Xq E Ai, there exist an open neighborhood W 
of Xq in M. and a continuous function a : W x W — ^ K such that (1) Cx{cr{x,y)) = y for all 
x,y eW (i.e. a{x,y) is a control taking x to y) and (2) for all x G W, cr(x,x) is the zero 
time control with T = and £i = Q for all i. 

Proof. Given a sequence of generators Yi, . . . , 1^, define the a sequence of compositions of 
formal exponentials by 



Q^iY^, ...,¥,):= e^^ o Q--\Y^, Y,^,) o e"^" o {Q''-\Y,, n_i))"\ (4c) 

If {Yi, . . . ,Y,y} are sufficiently close to zero, then using the Baker-Campbell-Hausdorff 
(BCH) formula [ ], it is easy to see that Q^^y^^y^) = el^^.m+jterms of order>2}_ 
Q^-i(yi,...,y;_i) = e^'^^-i -^dyayi+itermsoforderx.-i}^ ^-^^^ ^j^g formula implies that 

Q^iYi, ...,Y^)= e'^dy....ady2yi+{terms of order>i.}^ gy inductlou, this is true for all 1/ = 1, 2, ... . 

Now, define a corresponding sequence of maps -R^(lo) ^i? • • • ? ^u, ^1, ■ ■ ■ ,^u), where {^j} e 
M, by taking the sequence of exponentials in Q^(Yi, . . . ,Y^) {after having carried out the 
recursion in the above definition) and replacing e^^* by e^* ""i^oiCfcyfc fQj^ gome positive integer 
Then can be written, using BCH, as the exponential of a series of terms drawn 
from {^i^^^o; • • • ) ^^""^0; ^1^1; • • • ; ^vYy] and the Lie brackets thereof. Among the terms in 
this sum that do not include Yq, the lowest degree in the ^^'s is ^1 ■ ■ ■ adF^, . . .adF2^i, 
while the lowest degree term that includes Yq is degree 2ay in the ^^'s. So if we choose 

= \y/2 + IJ > Z//2, then 



Q\yi,y2) 



Y', y, -Y'2 -y 
e o e o e o e 



(4a) 
(4b) 



i?'^(yo,vi,...,n,ei,...,e.) = e' 



ad y,x-.. ad y2yi+{terms of degree in {^i}} 



(5) 
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where degree refers to the polynomial degree in the variables {^i, . . . ,Cu}, rather than the 
number of Lie brackets or vector fields involved. 

Using the fact that the Lie algebra generated by the control vector fields {fj}^i from 
(1) spans the entire tangent space Tx-M for each x, fix {X^j} C span{/i, . . . , /m} such 
that {Sfc := adXfci,^ . . . adXfc2Xfci} for k = l,...,n forms a basis for Txq-M. and define 
Fi : W X M'^ ^ by 



F,{x,r):= (/o,Xfcl,Xfc2,...,Xfc,,,sgnrfc|rfc|l/^^|rfc|l/'^^...,|rfc|l/'^^)j (x) (6a) 

gr„Sn+{higher order terms in r„} griHi+{highcr order terms in ''i}^^;) (6^) 

where the "higher order terms in r^" each have dependence of either \ri^\^^^^'^'' or r^lr^l^^'^'' 
for some positive integer j, so the higher order terms are functions of with deriva- 
tive zero at = 0. Then Fi is C^, 0) = x for all x G Ai, and |^(xo,0) = Hfc(xo) = 
adXkui, ■ ■ ■ adXfc2Xfei(a;o). Because these nested Lie brackets were chosen to be linearly inde- 
pendent, the Jacobian operator ^(xo,0) : M" — > T^o-M is full rank, hence an isomorphism. 
So by the implicit function theorem, there exists a neighborhood W C of a;o and a unique 
differentiable (i.e., C^) map : W x W — ?■ M" such that ip{x,x) = and Fi{x,(p{x,y)) = y 
for all X, ?/ G W. 

Now, observe that the sequence of exponentials in R'^(Yq,Yi, . . . ,Y^,C,i, . . . ,^^) can be 
interpreted as the solution of the control system (1) with piecewise constant controls. For 
example, in the case of -R^(/o, Yi, 1^2, ■Ci; ^2)5 where Yj = Yll^=i^jkfk^ the control C G K has 
final time T = + and control functions 



These controls can be shown to vary continuously with {^k} in the metric established for the 
control space IK (see Appendix C). Now define F : — )■ K to be the map tx{^F{r)) = Fi(x, r) 
that maps r G to the piecewise constant control in K whose integral curve corresponds 
exactly to the sequence of exponentials in the definition of Fi{x,r). Then F is a continuous 
map, and a : W x W — t- K defined by cr{x,y) = F{ip{x,y)) is also a continuous map. 
In addition, since ip{x,x) = 0, a{x,x) is the zero time control, because all of the ^^'s are 
zero, meaning that the final time T is zero. The map a is therefore exactly the desired 
cross-section map. □ 

We now return to the definition of a local fibration and define a function Aw : Xw — )■ PK, 
where Xyv C c^^(>V) x PW is the set of all (C,7) such that Cb(C) = 7(0), by letting 
A>v(C,7)(s) := C-ka{cb{C),'j{s)). Since Aw is continuous (see Appendix D) and (efe|E-i(Xw))° 
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Aw is the identity on Xw, Aw is a lifting function for ^b\^-i(^y^y Therefore tb is a local 
Hurewicz fibration and, by the Hurewicz uniformization theorem, also a global fibration. 
Note that Aw has the additional property that for a concatenation of unitary paths 71 ^72, 
Aw(C, 71*72) = Aw(Aw(C,7i),72)- 

B. Topology of the Fiber 

In order to establish the homotopy equivalence promised by Theorem 1, the following 
lemma is needed. 

Lemma 2. For any fibration p : E ^ B with contractible total space E and path con- 
nected base B, the map rj : VLB — t- F = p^^ipo) given by rj{^) = A(eo,7)(l) is a homotopy 
equivalence. Moreover, any map F — )■ QB that extends to a map g : E PB such that 
g{e){l) = p{e) is a homotopy inverse of rj, hence also a homotopy equivalence. 

Proof. Consider the homotopy H : PB x [0, 1] ^ E given by Hsi^j) = A(eo,7)(s), where cq 
is the basepoint of E and A is the lifting function of p. Let p : PB B he the endpoint 
map on the path space PB, so that p is also a fibration. Then Hi satisfies p o Hi = p, so 
that Hi is a fiber-preserving map, and Hq is the constant map -^0(7) = ^0 for 7 ^ PB. 
Let g : E ^ PB be any fiber-preserving map. Then g o Hi c:^ g o Hq which is a constant 
map on PB, and Hi o g c:^ Hq o g which is a constant map on E. Since PB and E are both 
contractible by assumption, go Hi and Hiog are homotopic to the identity maps on PB and 
E, respectively, so that Hi and g are homotopy equivalences. Moreover, since they are both 
fiber-preserving, they are fiber homotopy equivalences [27, §4H, Ex. 3] [2'.), §7.5], and fiber 
homotopy inverses of one another. Letting F denote the canonical fiber F := p^^(6o) and 
noting that the corresponding fiber of p is QB = p^^ipo), this implies that the restrictions of 
Hi to VLB and g to F are homotopy equivalences between these two spaces, and homotopy 
inverses of one another. □ 

In the case of the endpoint map of the control system C;, : K — )■ A^, for each s G [0, 1], 
let : K — i- K be the map Ps(T, £1, . . . , £m) = ((1 — s)T, Si, ... , Sm), where Sj(t) = £j{t) 
on [0, (1 — s)T] and is zero after (1 — s)T. Then ps{C) is continuous in C and s, po acts as 
the identity map on K, ps(0, 0, . . . , 0) = (0, 0, . . . , 0), and pi is the trivial map sending every 
control in K to the zero time control (0, 0, . . . , 0). So ps is a deformation retract of K to the 
point (0, 0, . . . , 0). Therefore K is contractible. Since Ai is assumed to be connected. Lemma 
2 implies that r] : VLAi — > F := e^^(6o) is a homotopy equivalence, where 77(7) = A(Co,7)(l), 
Co is the zero time control, and F C K is the set of controls taking the basepoint bo ^ Ai 
back to bo. Let r : K — )■ PB be the trajectory map, defined implicitly through the control 
system (1) with initial condition bo, that integrates the system to obtain the trajectory and 
then linearly reparametrizes the trajectory to obtain a curve in PB. Then p o t = p, so 
that r is fiber-preserving, so that, again by Lemma 2, r restricted to F is a homotopy 
equivalence between F and VLJ^. In addition, all fibers of are homotopy equivalent [27] 
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and the base point 60 is arbitrary, so Theorem 1 is proved. Furthermore, since is a 
second countable, hence separable, manifold, both Ai and QAi have the homotopy type of 
countable CW-complexes [oO] and therefore the fibers e^^(?/) all have the homotopy type of 
countable CW-complexes. 

Since e?, : K — )■ is a Hurewicz fibration with canonical fiber F := c^^(6), it admits a 
fibration sequence 

— > n'^K n^M ^riF ^riK^ riM 4 f A k ^ (s) 

where each triple X — )■ F — > Z is a fibration Y Z with fiber X up to homotopy equivalence 
[27]. The map i : F K. above is the inclusion map, and the map 1] : QAi — ?■ F is defined as 
in Lemma 2 [27]. For any pointed topological space {Z, zq), this fibration sequence gives rise 
to a long exact sequence of spaces of basepoint-preserving homotopy classes of basepoint- 
preserving maps [27] 

(z, n^K) {z, n^M) {z, nF) {z, nK) {z, nM) A {z, f) h {z, k) A (z, m). 

(9) 

The contractibility of K implies that fiK is also contractible, so that {Z, f2K) = {Z, K) = 
for all {Z, zq). Then the long exact sequence (9) contains the short exact sequence 

-> (Z, nM) ^ {Z, F) ^ 0. (10) 

Since control concatenation is continuous and the zero time control acts as the identity 
with respect to concatenation (i.e. 0-kC = C -kO = C), Fisan H-space [27]. This structure 
on F induces a natural monoid structure on {Z,F). In addition, the natural H-group 
structure on flAi induces a group structure on {Z,QM.). Since, as we noted above, the 
lifting function respects path concatenation, so too will r], i.e. r]{ui -k 0J2) = vi^i) ^^(1^2) 
for two loops uji,ijj2 G QJ^, and therefore 17* will respect the natural products on {Z,QJ\4) 
and {Z,F). Also note that the trivial loop (constant loop at I) in QA4 maps to the zero 
time control in F. Then the fact that 77* is a bijection (due to the exactness of (10)) means 
that the induced monoid structure on {Z, F) is, in fact, a group structure and rj^: is a group 
isomorphism. 

An obvious consequence is that all of the standard homotopy groups are isomorphic for 
these spaces. Indeed, iTiic'^iy)) := {S\c~^{y)) ~ {S\QM) =: 7ii{QM) ~ TTi+i{M) for each 
i > 0. In addition to this homotopy structure, the homotopy equivalence between c~^{y) 
and QAi implies that these two spaces share isomorphic homology groups [Hn{c~^{y)] G) ~ 
Hni^M] G)] and cohomology groups [if"(e~^(?/); G) ~ if"(f2A^; G)] for all coefficient groups 
G [27, Prop. 4.21]. We will take a closer look at these topological indices for the common 
cases of quantum control in the next section. 
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III. FIBER TOPOLOGIES IN QUANTUM CONTROL 



We now turn our attention back to quantum control and consider the detailed topologies 
of the fibers in the most common formulations, where either the manifold A4 is the special 
unitary group SU(A^) representing the space of unitary time propagation operators that 
are the general solutions of the Schrodinger equation, or the manifold M. is the complex 
projective space CP^^^ of states (modulo phase) of an A^-level quantum system, or the 
manifold is a generalized complex fiag manifold Fl(ni, . . . , n^) of unitarily accessible 
density matrices. As shown in Theorem 1, this reduces to understanding the topology of 
the loopspaces nS\J{N), nCF^-\ and ^]Fl(nl, . . . ,n^). 

A. Special Unitary Group 

First consider the case where the control problem ( 1 ) is a right-invariant affine system on 
SU(A^), i.e. the Schrodinger equation 

d 

ih-U{t)=H{t)U{t) t/(0)=I H{t)=H^ + Y,£iit)Hi (11) 

1=1 

where the ifj's are Hermitian and traceless and the control terms {iHi, . . . ,iHm} generate 
the Lie algebra su(iV). Theorem 1 states that for any target unitary propagator W G SU(A^), 
the set of controls in IK that take the identity I to W, c^-^iW), is homotopy equivalent to 
f2SU(A^). Now, the homotopy groups of SU(A^) can be shown by Bott periodicity [-31, 32] to 
be TTiiM) = for < z < 2iV- 1 and i even, 7Ti{M) = 0, and tt^A^) ~ Z for 2 < z < 2iV- 1 
and 2 odd. Therefore, 

i = 

1 < i < 2iV - 3 and i odd (12) 

2 < z < 2A^ - 2 and i even. 

As an immediate consequence, eif^(Vr) is both connected and simply connected for all N >2. 

The integral cohomology of i7SU(A^) was also investigated by Bott in the 1950's [32-34] 
and shown to be H*{QS'U{N); Z) Z[x2,a;4, . . . , X2N-2], the polynomial ring over Z, where 
the Xk are generators of degree k. Importantly, this implies that the cohomology of i7SU(A^) 
is trivial at odd dimensions. 



B. Complex Projective Space 

The space of pure quantum states of an A^-level system may be modeled as the complex 
projective space CP^~^, described by the fiber bundle — )• 5*^^"^ — )■ CP^~^, which gives 
rise to the long exact sequence 



vr,: c. 



9 



> 7r2{S') 712(5^^-') 7r2(CP^-^) ^ ni{S') ni{S^^-') ^ 7ri(CP^-^) ^ 0. 

(13) 

Since n^{S^) = for all i ^ 0, this implies that n,{CF^-'^) ~ ni{S^^-'^) for all i > 2. Also, 
since TTi{S^^'^) = for all i < 2A^ - 1, we have that Tri{CF^~^) = and 7r2(CP^-i) ~ Z. 



Therefore we find that for any two states ip, 
of t^^{(j)) are 



G CP , the first 2N — 2 homotopy groups 



VTi C, 





z 



z 



i = 

1 = 1 

2 < z < 2iV - 3 
i = 2A^-2. 



(14) 



The integral cohomology of fiCP^ ^ was also investigated in the 1950's [35, 36] and 
shown to be if*(fiCP^-^ Z) ~ /\[xi] (g) Z[?/i, ?/2, 2/3, •• • = (»i + j)!2/i+i} where Xi is 

a generator of degree 1 and yj is a generator of degree j * (2n — 2) for each j = 1, 2, . . . . 



C. Generalized Complex Flag Manifolds 

The closed-system evolution of a density matrix p (a positive, trace one, Hermitian ma- 
trix), can be described by the von Neumann equation 

, m 

ihj^p{t) = [H{t),p{t)] p(0) = Po H{t) = Ho + Y,£,{t)H, (15) 

i=l 

or, equivalently, p{t) = U(t)pQW{t), where U{t) G SU(A^) is the unitary propagator de- 
scribed in (11). Given some fixed initial condition po, the space of density matrices reachable 
by closed-system dynamics [assuming {iHj} generates su(A^)] is the orbit {UpoW : U G 
SU(A^)}, which can be regarded as the generalized complex flag manifold Fl(ni, . . . ,««) ~ 
SU(A^)/S(U(ni) © • ■ ■©U(nK)), where the n/s are the multiplicities of the eigenvalues of po- 
In the case where po is a pure state (i.e., po = |^o)(^o|)5 this is just the complex projective 
space, i.e. Fl(l,iV — 1) ~ CP'^"^. At the other extreme, if po is completely non-degenerate 
(representing a Boltzmann distribution over non-degenerate energy levels, for example), then 
the state space is the complete flag manifold Fl(l, 1, . . . , 1). 

These flag manifolds are compact, algebraic, homogeneous Kahler manifolds [ r"], si- 
multaneously possessing natural structures of complex analytic manifolds with Hermitian 
metrics, Riemannian manifolds, symplectic manifolds, and algebraic varieties. Topologically, 
they share the properties of being connected and simply connected, i.e. 7ro(Fl(ni, . . . , n^)) — 
TTi (Fl(ni, . . . , n^)) — 0. Also, from the homotopy long exact sequence of the fibration 

S(U(ni) © • • ■ © U(n,)) ^ SV{N) ^ Fl(ni, . . . ,n,), (16) 
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it may be seen that 7i2{F\{ni,. . . ,n^)) ~ 7r2(S(U(ni) © ■■■ ©UK)) ~ Z^-^ Therefore 
the loop spaces are all connected, 7ro(fiFl(ni, . . . ,7^^)) — 7ri(Fl(ni, . . . jn^)) — 0, and have 
nontrivial fundamental group: 7ri(f2Fl(ni, . . . , n^)) — 7r2(Fl(ni, . . . , n^)) — Z**^^. Higher 
order homotopy information about these spaces will depend on the choices of the parameters 
{rii, . . . , n^), i.e. on the eigenstructure of po- 

The integral homology ring of the loop space of the complete flag manifold was just 
recently shown to be 

if*(fiFl(l,l,...,l);Z) ~ H,{Q{S\J{N)/T^-^);Z) 

~ (T(xi,. . . ,XAr__i) © Z[?/i, . . .,yN^i])/{xl = XpXg = 2yi for 1 < k,p,q < N -1, p ^ q), 

(17) 

where the generators xj are of degree 1, and the generators yj are of degree 2j, and 
T(xi, . . . , xn-i) is the tensor algebra generated by xi, . . . , xn-i [38]. The integral (co)homologies 
of the loop spaces of the other generalized complex flag manifolds appear not to have been 
published, but both the rational and mod p cohomologies have been investigated [39]. 



IV. TOPOLOGY OF LANDSCAPE LEVEL SETS AND CRITICAL SETS IN 
QUANTUM CONTROL 

In the study of quantum control landscapes, we are often interested in "dynamical" 
subsets c'Ij.^{A) C K corresponding to "kinematic" subsets A d M.. Most typically, these 
kinematic subsets are level sets and critical sets of objective functions defined on A^. In this 
section, we consider the tools that may be used to investigate the topology of these types 
of dynamical subsets and analyze a few examples from quantum control. In particular, we 
will examine the topologies of level sets and critical sets of two classes of quantum control 
landscapes. 

The first of these families of landscapes consists of the unitary quantum control land- 
scapes J{U^ = I Tr(ylA"''l^"'"[/)p where A is any complex N x N matrix and is a target 
unitary transformation (i.e., quantum gate). These landscapes generalize the standard fi- 
delity measure used in quantum information applications. The critical point set of such a 
landscape in the case where A is non-singular can be shown (aside from the J = min- 
imum set) to comprise disjoint manifolds, each of which is a direct product of complex 
Grassmannians [10, 40]. If AA^ has eigenvalues with multiplicities rii, . . . ,11^, then these 
critical submanifolds are of the form Gr^, (C"i )©■■■© Gr^^ (C"'=) for some set of parameters 
< i^j < Uj. It should be noted that for any such critical submanifold C C SU(A^) and any 
root of unity e^'^*^/^, the set e^'^'^^/^C is another critical submanifold, disjoint from C, with 
the same critical value as C. 

The second landscape family can either be described on the special unitary group by 
J{U) = Tt{UpoWO) or on the unitary adjoint orbit through po [i-e. the flag manifold 
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Fl(ni, . . . , n^)] by J'{p) = Tr(p(9). In either case, O is a Hermitian matrix describing 
a particular quantum mechanical observable, and the landscape represents the value that 
would be measured on a quantum system with state UpoU^ (or p, respectively). It may be 
shown that the natural action of SU(A^) on F1(?t,i, . . . , n^) is a Hamiltonian action and that J' 
is the associated "Hamiltonian function" that generates this action and defines the moment 
map. As such, J' is a Morse-Bott function such that each level set is connected, and therefore 
J' has exactly one local minimum submanifold and one local maximum submanifold, i.e. J' 
has no "traps" [41, 42] [43, §3.3-3.4] [44, Ex. 3.24]. Furthermore, the critical submanifolds 
of J' are diffeomorphic to direct products of generalized complex flag manifolds: C = 
Fl(/cii, . . . , kis) © ■ ■ ■ © F1(A;ki, • • • , ^ks) where J2'j=i ^ij = ''^i each i. As is the case with 
the adjoint orbit itself, these critical submanifolds are homogeneous Kahler manifolds [37]. 
When formulated on the unitary group U(A^), J was shown to have critical submanifolds 
diffeomorphic to (U(n) ©U(m)) /Gtt, where is a subgroup of the form {V (BttVtc^ : V G 
U(n) and tiVti'^ G U(m)} for some permutation matrix vr, U(n) = U(ni) © ■ ■ ■ © U(nK) and 
similarly for U(m), and where n = [ni, . . . , ««]""" and m = [mi, . . . , TUrY are the multiplicities 
of the eigenvalues of p and C, respectively [ ]. This is the fiber bundle defined by the 
restriction of / : U(A^) — Fl(r;,i, . . . , ra^) to /^^(C) with base C a critical submanifold 
of J' and fiber U(?2i) © ■ ■ ■ © \]{n^). Likewise the critical set of J over SU(A^) comprises 
submanifolds diffeomorphic to the fiber bundles given by the restriction of g : SU(A^) — >■ 
Fl(ni, . . . ,n«;) to g'^{C') with base C and fiber S(U(ni) © • ■ ■ ©U(n«;)). 



A. Homotopy 

For any subset A oi M.^ restricted to t~^{A) is a Hurewicz fibration with the same 
fiber as (up to homotopy equivalence) [27]. So assuming that A is connected, one way to 
examine the topology of t~^{A) is through the long exact sequence 

■■■^7r2(e;^(yi))^7r2(/l)^7ri(F)^7ri(e;i(^))-^7ri(yi)^7ro(F)^7ro(e;i(yi))->0, (18) 

which, because the fiber F = c^^{x) (the set of controls steering from x back to x) is 
homotopy equivalent to QAi by Theorem 1, may be rewritten 

>7r2(e;^(yi))^7r2(yi)^7r2(A^)^7ri(e;^(yi))^7ri(yi)^7ri(A^)^7ro(e;^(yi))^0. 

(19) 

This long exact sequence may be used to obtain higher order homotopy information about 
c~^{A). For example, when the fiber F is connected and simply connected (as is the case 
when Ai = SU(A^)), the fundamental group 7ii{t~^{A)) is isomorphic to 7ri(>A). On the 
other hand, when Ai = CP^~^, the many trivial homotopy groups of the fiber indicated in 
(14) imply that TTi{c~^{A)) ~ TTi{A) for alH = 3, . . . , 2A^ — 3. Much information about the 
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structure of these dynamical sets c~^{A) can be obtained through a careful analysis of this 
long exact sequence and utilizing other tools of homotopy theory, including perhaps spectral 
sequences. The implications for connectedness of these sets will be considered presently, 
however an exhaustive analysis of the higher homotopy structure is outside the scope of this 
paper. 



B. Connectedness 

An immediate consequence of the homotopy long exact sequence (18) is that if A is 
connected and ni^A) — )■ 7ro(F) is surjective, then c~^{A) is connected. This map is defined 
in a similar fashion to t] from (8): given an element of 7ri(yi), choose a representative 
loop in A, lift it to t~^{A) by the lifting function, and select the connected component of 
the endpoint of the resulting curve to obtain the corresponding element of ttq{F). Since 
7ro(F) ^ ttq^QJH) ~ 7ri(A^), surjectivity of the map 7ri(yi) — )• 7ro(F) is equivalent to the 
statement that every element of the fundamental group of A4 has a representative loop in 
A, or that : 7ri(yi) — t- 7ri(A^) is surjective, where i : A ^ Ai is the inclusion map. More 
generally, if tti^Ac) — )■ tto{F) is surjective for each connected component of Ac C A, then 
c~^{A) has exactly the same number of connected components as A. In particular, if F is 
connected then this is satisfied trivially. Since the fibers of the fibrations discussed in the 
last section are all connected, this statement is especially relevant to quantum control. 

Since each level set of J'{p) = Tr(p(9) is connected in Fl(ni, . . . , n^), and Fl(ni, . . . , n^) 
is simply connected, the above arguments imply that the corresponding "dynamical" level 
sets in IK are also connected. Of course, this conclusion is independent of the choice of state 
space. In other words, had the control problem been defined on = SU(A^) with the 
kinematic landscape J{U) = Tt{UpoWO), the dynamical level sets are identical to those of 
the Fl(ni, . . . ,72^) construction, and therefore are also connected. 

With the unitary landscapes J{U) = \ Tr{AA^W^U)\'^ on SU(A^) where A e C^^^ is 
an arbitrary parameter matrix [10, 10], the connectedness story is different. When A is 
nonsingular, the maximal level set of J comprises exactly N discrete points, one for each 
global phase rotation of by an A^'th root of unity, i.e. A = {e^^W : /c = 0, . . . , A^— l} C 
SU(A^). As a result, every level set from that maximal value down to the next critical value 
consists of exactly A^ connected components (homeomorphic to spheres). And therefore, each 
of the corresponding "dynamical" level sets comprises A^ connected components. This result 
is perhaps better understood in terms of control on the projective unitary group PU(A^). 
Since the global phase of the unitary evolution operator (or indeed the quantum state) is 
not physically measurable, the space of evolution operators should really be considered to be 
U(A^)/U(1)I = PU(A^). This is a Lie group with fundamental group 7ri(PU(A^)) = Z/NZ, 
i.e. the additive group of integers mod A^. Since the landscape J is invariant to global 
phase, it is well-defined on PU(A^) and the maximal level set of controls is the same as 
when the control problem is defined on SU(A^). But the interpretation now is that, because 
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7ri(PU(A^)) = Z/A^Z, for any given target W there are exactly homotopy classes of 
paths joining the identity I to the target W, and each class is uniquely identified with one 
connected component of the maximal level set of J in K. This then leads to the possibility 
of a physical interpretation of control set topology because it may be reasonable to identify 
these N homotopy classes of paths with N classes of physical mechanisms of control for 
arriving at the target operator W. Then each connected component of the maximal level 
set in K would uniquely correspond to one of these N mechanism classes. 



C. (Co)Homology 

Since restricted to c~^{A) is a Hurewicz fibration with the same fiber as c^, the most 
promising method for evaluating the (co)homology of the dynamical set e|j^^(7l) is the Serre 
spectral sequence. The Serre spectral sequence applies to Hurewicz fibrations and attempts 
to compute the (co)homology of the total space [in this case e~^(yi)] given the (co)homologies 
of the fiber F ~ QAi and the base A. The sequence appears in various forms, the simplest of 
which is dependent on the satisfaction of one additional assumption, namely that the action 
of the fundamental group of the base A on the (co)homology of the fiber is trivial [45, 46]. 
Since the state spaces we have considered - SU(A^), CP^^^, and Fl(ni, . . . ,71^.) - are all 
simply connected, every loop in these spaces acts trivially on the (co) homology of the fiber. 
For any subset A of one of these spaces, a loop in yi is a loop in the state space and therefore 
acts trivially on the (co)homology of the fiber. So, even though A may not be simply 
connected, its fundamental group acts trivially on the (co) homology of the fiber simply by 
virtue of A being a subset of a simply connected space, and the fibration being defined 
by restriction to t~^{A). As a result, the Serre spectral sequence beginning with ^ 
Hp{A;Hg{F;G)) [or E^^ = Hp{A; Hi{F;G))] converges to H,{t-\Ay,G) [H* {c-\Ay, G) , 
respectively] for every abelian group G of coefficients [46, Thm. 5.4] [45]. 

While this gives the starting point of the spectral sequence (assuming the (co)homology of 
A is known), the effort involved in working out the limit of the sequence can be substantial. 
For this reason, a full accounting of the (co)homologies of, say, the dynamical critical point 
sets of Tr(pC), is beyond the scope of this paper. However, in certain special cases, the 
Serre spectral sequence converges immediately at the E"^ step. In particular, this occurs 
when the (co)homologies of both the fiber F and the base A are torsion free and zero for 
odd dimensions. Under these conditions, the integral homology of the total space c~^{A) 
is simply given by H^{t~\A)]Z) = H^{A;Z) ®if,(F;Z), and similarly for the integral 
cohomology [ i5]. 

This enables a vast simplification of the problem of computing the integral cohomol- 
ogy of the dynamical critical sets of | Tr(y474^iy^?7)|. The integral cohomologies of the 
complex Grassmannian manifolds Grj,(C") and of the loop space of the special unitary 
group f2SU(A^) are free of torsion and zero for odd dimension [ ] as desired for the sim- 
plification indicated above. Because a critical submanifold C of | Tr(yly4"l'iy"''[/)| in SU(A^) 
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is a direct product of Grassmannians, the Kiinneth formula may be invoked to conclude 
H*{C; Z) = H* {Gi^, (C"i)) (g) ■ ■ ■ (g) if* (Gr^^ (C"«)) . It follows that the integral cohomology 
of a dynamical critical set corresponding to the preimage of a critical submanifold C of 
I Tr(y4y4^iy^f/)| in SU(A^) is torsion-free, zero for odd dimensions, and is given by 



H*{ci\Cy,Z) =ii*(Gr,,(C"i)) ®---®if*(Gr,^(C"-)) ®Z[x2,...,X2iv-2], (20) 

where the X2j are free generators of degree 2j, and the cohomologies of the individual 
Grassmannian manifolds are given by the Poincare polynomial [ 17] 

Pi[Gr4C")] = ^^_^2)...(l_i2.)(l_i2)...(l_^2(n-.))- (21) 



V. CONCLUSIONS 



A proof of a generalization of Sarychev's theorem was presented, demonstrating that the 
space of controls driving the dynamics of an affine control system from a given starting point 
to a given target point is homotopy equivalent to the loop space of the state manifold of 
the system, provided the non-drift part of the control system satisfies the Lie algebra rank 
condition. Moreover, the proof shows that the endpoint map is a Hurewicz fibration. 
These results enable the analysis of the topological structure of subsets of control space 
of the form z^^{A) C K. Both the homotopy and (co) homology structures of critical sets 
and level sets may be examined using tools of algebraic topology such as exact sequences 
and spectral sequences, as demonstrated in the simplest cases. It should be stressed that 
these results depend heavily on the control space being considered. In particular, the proof 
does not apply in the case of scalar control and/or a fixed final time T, which have often 
been the conditions in prior quantum control landscape analyses. However, in the case of 
a two-level quantum system, Montgomery [22] has essentially shown how Smale's work on 
regular homotopy [2 1] can be applied to a scalar control problem on S0(3) ~ PU(2), leading 
to a similar conclusion to Theorem 1 (see Appendix E). Similar analysis for scalar control 
problems on higher-dimensional quantum systems may be more difficult, possibly due in 
part to the presence of singular controls [48]. 

The most straight-forward topological information that can be obtained is the connect- 
edness of these sets. It was shown that, since the state spaces of greatest interest in closed 
system quantum control are generally simply connected, a "dynamical set" of the form 
t~^{A) has exactly one connected component for each connected component of the "kine- 
matic set" A. This implies, for example, that the dynamical level sets of the ensemble 
landscape Tr(p(9) corresponding to a quantum mechanical observable O are all connected 
and that the highest dynamical level sets (from the maximum value down to the next highest 
critical value) each comprise exactly N connected components for an A^-level system. This 
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structure of the dynamical level sets and critical sets can be important for understanding the 
behavior of optimal control algorithms employed both in the laboratory and in numerical 
simulations, especially continuous level set and critical set exploration algorithms such as 
D-MORPH [11-15]. In some cases, for example the case of the set of all controls producing 
a target unitary evolution operator up to global phase, the connected components of the 
control set may also have a physical interpretation in terms of control mechanism or other 
characteristics. In addition, this connectedness story as well as the wealth of additional 
homotopy and (co)homology content of the theorem offer a tantalizing glimpse of a deeper 
structure of dynamical quantum control landscapes and more generally of quantum optimal 
control theory, the full implications of which remain to be understood. 
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Appendix A: Continuity of the Endpoint Map 

In this appendix, we prove the claim that the endpoint map e^^ : K — )■ is continuous. 
From [_ , Theorem 10, Chapter 4] comes the result that for an affine control system on a 
Riemannian manifold, convergence of controls in L^[0,T] (for a fixed final time T) implies 
uniform converges of the trajectories and therefore convergence of the endpoints. Suppose 
that, in K, Ck = {Tk, . . . ,8^) -> {T,Si, . . . , Sm) = C. Then ^ T, and in ^^[0, f ] 
for T = sup{Tfc} < oo, we have that £^ — )■ £j for each j. Then the trajectories converge 
uniformly Zk{-) — )■ z{-) on [0,T]. Each Zk is the controlled trajectory over [0,Tfc] followed 
by pure drift (corresponding to zero control) on [Tfc,T]. Then supjg (i(2;fc(t), z(t)) — >• as 
A; — )■ oo, where d{-,-) denotes the topological metric on M. induced the Riemannian metric. 
Now, 

sup d{zk{TTk),z{TT)) < sup d{zk{TTk),z{TTk)) + sup d{z{TTk), z{tT)) (Ala) 

re[0,l] -re[0,l] -rG[0,l] 

< sup d{zk(t), z(t)) + sup d{z{TTk), z{tT)). (Alb) 

te[o,f] re[o,i] 

The first term of (Alb) vanishes as — oo due to uniform convergence of to z on [0, T]. 
Since z is continuous on [0, T], by the Heine (-Cantor) theorem it is uniformly continuous, so 
that the second term also vanishes as /c — )■ oo and — ?■ T. As a result, the trajectories Zk on 
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[0,Ta;], linearly reparametrized to [0,1], converge uniformly to 2;, similarly reparametrized. 
This means first that the map from K to this linearly reparametrized trajectory space with 
uniform topology is a continuous map, so the metric topology we have defined on K is strictly 
stronger than this uniform topology on trajectory space. And second, that the endpoint map 
is continuous. 



Appendix B: Continuity of Control Concatenation 

We defined concatenation of controls in (3) with notation C -k D ioT C,D G K. In this 
section we prove that this operation 7ir:KxK— i-Kisa continuous map. Since K is a 
metric space, K x IK is also metric, so it suffices to consider the convergence of Ck-^Dk when 
{Ck,Dk) — )■ {C,D) and therefore ^ C and Dj. — ?■ D. We will write the components of 
these controls as C = (T^, £1, . . . , £m), Ck = (Tf , ...,£lj), D = (T^, J^i, . . . , J^m) and 
Dk = {T^, J^i, . . . , J^m)- By the triangle inequality, d{Ck * -Dfc, C ^ D) < d{Ck * Dk, Ck * 
D) + d{Ck -k D,C -k D). It may easily be observed from the definition of -k and of the metric 
d on K, that d{Ck ^ Dk, Cki< D) = d{Dk, D). Now 



f|Tf-T^| + E,^i/o^'l^^W-^.W|dt 



d{Ck*D,Ci<D) = < 



+ \£;it)-J^,it-T^)\dt 
+ /;°o|^,(t-Tf)-J-,(t-T^)|dt T^<T, 



(Bla) 



+ /^o"|^,(t)--F,(t-Tf)|dt 

+ /;^e|j-,(t-Tf)-j-,(t-r^)|dt Tf<r^ 



<d{Ck,c) + y2 / iJ'Mdt 

+ / |j;.(t + |Tf -r^|)-j-,(t)|dt. 
^0 



(Bib) 



Since the continuous functions are dense in L^, let {Qj} be continuous functions supported 
on [0,T^ + 1/^] such that Q'- Tj on [0, 00). Then for each /, 



lim / |JS-(t)| dt < - Q\\\ + lim / \Qit)\dt = ^ - (B2) 

For any e > there exists L such that for all I > L, \\J^j — < e, so we see that the left 
hand side of (B2) is zero. Likewise, the third term of (Bib) can be bounded for each / as 



lim 



J-,(t+|Tf-T^'|)-J-,(t)|dt < 2||J-,-e;'l|+lim 



^j(t+|Tf-T^|)-^j(t)|dt. 

(B3) 
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Since Qj is a continuous real-valued function on a compact interval, by the Heine (-Cantor) 
theorem, it is uniformly continuous, so for each e > 0, there exists K such that for all k > K, 
\g^j{t -f |Tf - T^l) - < e for all t e [0, + 1//]. So the integral on the right hand 

side of (B3) is bounded by e(T^ + 1/0 arbitrarily small e > and so it vanishes in 

the limit as — )■ cxd, meaning that the left hand side is bounded by 2||J^,• — Qj\\ for any /. 
Since this can be made as small as desired, the left hand side of (B3) must vanish. So we 
finally conclude that 

d{Ck i<Dk,Ci<D)< d{Ck * Dk, Ck*D) + d{Ck * C ^ L>) ^ (B4) 

so that concatenation is a continuous operation in the K metric. 

Appendix C: Uniform Continuity of the Cross-Section Map 

Recall from Lemma 1 that the cross-section map (TiWxWcA^— t-K was defined 
to be a{x,y) = F{x,ip{x,y)), where W is an open neighborhood of xq G Ai, F : — )■ K, 
and (y9 : W X W — )■ M"'. Since (p was defined and shown to be by the implicit function 
theorem, the proof that a is continuous reduces to showing that F is continuous. Then 
uniform continuity follows from observing that, by restricting to a smaller neighborhood of 
Xq if necessary, a is continuous on W x W, which is compact. So invoking Heine (-Cantor), 
we can conclude that a is uniformly continuous. 

Turning now to the proof that F is continuous, let us first consider maps 

i?'^(/o,ri,...,n,6,.--,U = {T,s,,...,SM)eK (ci) 

that produce the piecewise constant controls whose integral curves are exactly the sequence 
of exponentials in the definition of W^, hence z o = W^. Then there exist 5*^ : — )■ K 
for = 1, . . . , 2-^ + 2''-^ - 2, such that for fixed Yi, . . . ,Y^, R'' = Si ^ S2 * ■ ■ ■ 82^+2-^^2, 
and where each Sk{^) = {T^, J-'f , . . . , J-'^j) depends on only one of the ^j's, is defined on an 
interval of length = ^j""" , and each J^f takes on a constant value proportional to ^j"^""". 
Now, it is easy to see that each Sk is continuous in ^. So, by the continuity of -k proved in 
Appendix B, i?^ depends continuously on ^. Then, since F{r) is formed by concatenating 
the output of these R'^ functions with the input of each R'^ depending continuously on some 
yj, we have that F is continuous. And by the above arguments, we conclude that a is 
uniformly continuous. 

Appendix D: Continuity of the Lifting Function 

In this appendix, we prove that the (local) lifting function Aw : Aw — ^ -PIK is a continuous 
map in the appropriate topologies. Recall that W C is an open set and Aw C c^^(>V) x 
PW consists of pairs (C, 7) such that Ch{C) = 7(0). PW is given the relative topology 
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as a subset of the path space PM., which itself is given the compact-open topology. Since 
is a metric space under the bi-invariant topological metric induced by the Riemannian 
metric from the real Hilbert-Schmidt inner product, both PAi and PW are also metric 
spaces. Likewise PK, c~^(>V), and Xyv are all metric spaces. So to prove continuity of Aw, 
it suffices to prove convergence of A>v(Cfc,7fc) to Avy(C, 7) when {Ck.lk) (C*)7)- 

Recall that A^ was defined by Aw(C, 7) = C*cr(ef,(C), 7(5)), where a is the cross-section 
map defined in Lemma 1 whose uniform continuity was established in Appendix C. By the 
triangle inequality and the definitions of the metric on K and of the concatenation operator 



ci(Aw(Cfc,7fe), Aw(C,7)) = sup d[Ck-k a{tb{Ck),lk{s)),C -k a{tb{C),-i{s))) (Dla) 

s6[0,l] 

< sup (i(Cfc^cr(e6(Cfc),7fe(-5)),C'fe*o-(efe(C),7(s))) 

sG[0,l] 

+ sup d{Ck k a(e,(C), 7(s)), C * a{z,{C),^{s))) (Dlb) 
se[o,i] 

= sup d{a{tb{Ck),lk{s)),ci{zb{C),-i{s))) 

se[o,i] 

+ sup d(Cfc^(T(e5(C),7(s)),C^a(efe(C),7(s))). (Die) 

se[o,i] 

Consider the first terms on the right hand side of (Die). Since by assumption 7^ — )■ 7 
uniformly and Ck C, it follows by continuity of Cb that Cb{Ck) — )■ Cb(C), and therefore 
(2fe(Cfc), 7fc(s)) — )■ (efc(C),7(s)) uniformly over s G [0,1]. Then the uniform continuity of a 
implies that o'{tb{Ck),'yk{s)) — t- cr(e5(C), 7(5)) uniformly, so this first term converges to zero 
as — i- 00. 

Turning now to the second term on the right hand side of (Die), observe that {Ck} is a 
compact subset of K since any open cover must admit an open set containing C and this set 
will also contain all but finitely many of the C^'s, so that only finitely many additional sets 
from the cover are needed to construct a finite subcover. Also {cr(e6(C), 7(5)) : s G [0, 1]} 
is a compact subset of IK since it is the continuous image of a compact set. Therefore 
A := {Ck} X {(T(efe(C), 7(5)) : s G [0,1]} is compact in K x K. Since -k is continuous 
over K X K, it is uniformly continuous over A by the Heine (-Cantor) theorem. So for any 
e > there exists a. 6 > such that for any (ai,/3i), (0:2, /32) £ A with d{ai,a2) < S and 
d{f3i, (32) < S, d{ai -k (3i,a2 * P2) < e)- Then for any e > 0, there exists K such that for 
all k> K, d{Ck,C) < 6, and therefore d(Cfe ^ (T(7(s)e(C)t), C * a(7(s)e(C)t)) < e for all 
s G [0, 1]. Therefore the second term on the right hand side of (Die) converges to zero, so 
the left hand side converges to zero, and therefore the lifting function Aw is continuous with 
respect to the given topologies. 
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Appendix E: Dynamic Homotopy of Scalar Control on 2-Level Systems 



Consider a right-invariant affine control system of the form 

i'^ = {H^ + S{t)H;)U{t) t/(0)=I (El) 

on the 2x2 special unitary group SU(2), hence Hq and Hi are Hermitian and trace zero. 
We will assume that the system satisfies the Lie algebra rank condition (LARC), so that 
[iHQ^iH^ lies outside the linear span of {iHq.iHi}. In this appendix, we describe the 
relationship between the topological study of the trajectory space of this scalar control 
problem and Smale's work on regular homotopy. 
For a given control let 

Ao := ^ A{t) := U\t)A^U{t) (E2a) 
ll-"i|| 

^° ■= WTr^\ ■= U\t)B,U{t) (E2b) 
\\[iHq,iHi\\\ 

^" - u\t)C,u(t) (E2c) 

\\iHq - {iHq,Aq)Aq\\ 

so that A{t), B{t), and C{t) are curves on the unit sphere within the three dimensional su(2). 
Moreover, Aq, Bq, and Cq may be seen to comprise an orthonormal basis for su(2) under 
the Hilbert-Schmidt inner product, so that A{t), B{t), and C{t) are a rotating orthonormal 
frame. Then it is readily found that 

A A 

— it) = U\t)[tHit),Ao]Uit) = aB{t) (E3a) 

^(t) = U\t)[iH{t),B^]U{t) = -aAit) + {l3 + ^£{t))C{t) (E3b) 
An 

— it) = U\t)[tH{t),Co]U{t) = -{p + ^S{t))B{t) (E3c) 

where 



/3 := {zHo,Ao)\\[Co,Ao]\\ 7 := ll[Co,^i^i]||. (E3d) 

So A{t) travels at constant speed a over S"^ C su(2) with continuous velocity vector aB(t). 
Therefore, A{t) is a regular curve [ ] on 5*^ C su(2). The covariant derivative of B{t) is 
just ^ = {f3 + j£{t))C{t) so that B{t) can change arbitrarily quickly in either direction (C 
or —C) by tuning the control S{t). As a result, any regular curve in S"^ C su(2) starting at 
Aq with initial normalized velocity Bo and with constant speed a may be generated in this 
way by choosing the correct control function S and final time T. In other words, since the 
curves A{t), B{t), and C{t) uniquely determine the propagator in PU(2), there is a one-to- 
one correspondence between trajectories of (El) on PU(2) and regular curves on S"^ G su(2) 
starting at Aq (up to reparametrization) . This construction is essentially the same as the 



a := \\[iHo,Ao] 
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example described by Montgomery [22] relating control trajectories on S0(3) (~ PU(2)) to 
regular curves on 5*^. 

Stephen Smale studied regular liomotopies of regular curves on Riemannian manifolds 
[21]. In particular, he looked at the set E{pi,vi,p2,V2) of regular curves starting at a point 
Pi with initial velocity vector vi and ending at a point p2 with final velocity V2- He derived 
a weak homotopy equivalence between E{pi,vi,p2,V2) on a manifold Ai and flTQ^Ai), the 
loop space on the unit tangent bundle of the manifold, where the unit tangent bundle 
To(A^) is the bundle of all unit length (i.e. normalized) tangent vectors to Ai. Now, the 
unit tangent bundle of the sphere, To(S'^), can be identified with S0(3) by thinking of 
the first column of O G S0(3) as the point in S"^, the second column as an arbitrary unit 
tangent vector to S"^ at that point, and the third column as simply the cross product of the 
first two, carrying no additional information. Thus, since S0(3) ^ PU(2), Smale's theorem 
shows that E{A{0), B{0), A{1), B(l)) (where curves have been reparametrized to travel at 
constant speed on [0, 1], eliminating the final time T as a variable) on S"^ C su(2) is weakly 
homotopy equivalent to f2PU(2), the loopspace on PU(2). 

It may be observed that the topology given by Smale to the space of regular curves on 
is identical to that induced by the compact-open topology on the path space of PU(2), 
so that the one-to-one correspondence established above between trajectories of the control 
problem on PU(2) and regular curves on S"^ is a homeomorphism between these spaces. This 
result implies the existence of a weak homotopy equivalence between the set of trajectories 
of (El) joining I e PU(2) to a target W E PU(2) and the loop space of PU(2). Thus, we 
have the following theorem: 

Theorem 2. Given any W G PU(2), the space of all trajectories from I to W of the control 
system (El) is weakly homotopy equivalent to QP\J{2), the loop space o/PU(2). This loop 
space (and therefore the indicated set of trajectories) comprises two connected components. 

Moreover, using the fact that SU(2) is the universal covering group for PU(2), it is 
straightforward to lift this result to SU(2), yielding the theorem 

Theorem 3. Given any W G SU(2), the space of all trajectories from I to W of the control 
system (El) is weakly homotopy equivalent to f2SU(2), the loop space o/SU(2). In particular, 
this loop space (and therefore the indicated set of trajectories) is connected. 

In addition, since Smale has shown [2 i] that the endpoint map on the space of regular 
curves on 5*^ is a Serre fibration [27] , and since this space of curves in homeomorphic to the 
space of PU(2) trajectories of (El), the fact that the projection vr : PU(2) — )• CP^ onto the 
Riemann sphere CP^ (i.e. the Bloch sphere; i.e. the complex projective line) is a fibration 
may be used to show that the endpoint map on the space of CP^ trajectories for some 
initial state po 7^ 1/2 is also a Serre fibration. Similar to the proof of theorem 1, using the 
contractibility of this space of trajectories, this implies 

Theorem 4. Given any 2x2 density matrix po 7^ 1/2, and any pi isoentropic to po (i.e. 
possessing the same eigenvalues), the space of all trajectories from po to pi of the control 
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system 

t^ = [Ho + £muPit)] p(0)=Po (E4) 

at 

is weakly homotopy equivalent to flCP^, the loop space ofCP^ ^ S^. In particular, this loop 
space (and therefore the indicated set of trajectories) is connected. 

It is speculated that these resuhs could be pulled back to control space as in Theorem 
1, however this would require redoing Smale's arguments from the control perspective and 
with a specific topology on the control space in mind. 
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